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Abstract 

We present an elementary algorithm for the dynamics of recently introduced soli- 
ton cellular automata associated with quantum afiine algebra Uq{g„) at g = 0. For 
0n = ^4^^'', the rule reproduces the ball-moving algorithm in Takahashi-Satsuma's 
box-ball system. For non-exceptional g„ other than A^^', it is described as a motion 
of particles and anti-particles which undergo pair-annihilation and creation through 
a neutral bound state. The algorithm is formulated without using representation 
theory nor crystal basis theory. 



1 Introduction 

The box-ball system invented by Takahashi and Satsuma Jt^ , ^ is a remarkable example 
of one dimensional soliton cellular automata. It represents the dynamics of balls hop- 
ping along the array of boxes under a certain rule. The evolution equations have been 



studied extensively in |TTMS, TTM, TNS| by means of the ultradiscretization of soliton 
equations. 



In llHKTll , |FOY| , |HIIIKTT[ , the box-ball system was identified with a solvable vertex 
model in statistical mechanics Q in a low temperature limit. Mathematically it brought 
a crystal basis theory Q into the game, and led to a generalization lab eled wi th afSne 
Lie algebras g„. We call the resulting sy stem the fln-automaton. See | HKT1 for the 
construction for non-exceptional g„, and [EKOTYj for the study of soliton scattering. 
The box-ball system corresponds to the g„ = An^ case | FOY , HIIIKTT |. 



The dynamics of the g„-automaton is governed by the combinatorial R |KMN] in the 
crystal theory, which specifies the local interaction of a box and a carrier in the language 
of the box-ball system. Except the A^P case |NY], the complexity of the combinatorial 



R HKOTl , HK0T2| makes it tedious to compute th e dynam ics of the associated g„- 
automaton. Fortunately the difhculty was overcome in |HKT2| for the automaton having 
a carrier with infinite capacity, where a factorization of the time evolution into simple 
Weyl group reflections was proved. 



The aim of this paper is to translate the factorized dynamics | HKT2 1 stated in the 
crystal language into an analogue of the ball-moving algorithm, and thereby to describe 
the g„-automaton in terms of a simple evolution rule. For simplicity we shall restrict 
ourselves to the case where all the boxes have minimal capacity. In the crystal formulation 
of the time evolution: Bm ® (■ • ■ ® B^i <8) -Bmi+i ® • ■ • ) ~ (■ • ■ (8) B^i ^ -Bm.+i ® ■ ■ • ) ® ^m, 
the situation corresponds to the choice Vmi = 1 and Af = oo. The factorized dynamics 
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for ^l^-* -automaton reproduces the decomposition of the ball- moving algorithm in the 
box-ball system into a finer procedure where one only touches the balls with a fixed 
color. For the other g„'s in question, the dynamics is decomposed similarly, where each 
finer procedure (denoted by Kj in Section ^) is now described as the motion of particles 
and anti-particles with a fixed color exhibiting pair-annihilation and creation through a 
neutral bound state. 

In Section || we specify the set of local states and give the algorithm. In Section ^ we 
present examples. In Section |^ we list the patterns that behave as solitons. In Section 
H we briefiy explain the relation of the algorithm in Section |^ and the result in [ HKT2{ . 
The key is the formula (^). 

Although our notations originate in the crystal theory, the contents in Sect ions 
use no result from it. Despite being just a translation of a more general result in [HKT 
we hope this paper possesses its own role to make the B„-automaton more accessible and 
familiar to the people working on cellular automata and discrete integrable systems. 



2 Algorithm 



We formulate the g„-automaton associated with the non-exceptional affine Lie algebra 
g„ = AiP{n > l),4l-i(" > 3),4«(n > 2),B^^\n > 3),dn\n > 2),Dk^\n > 4) 

(2) 

and > 2). Actually the algorithm itself given below is meaningful formally for 

A'^\ B2^\ 0^2"^ ^iid Dg^-* as well. First we specify the set B of local states and the 
sequence ji, . . . , jrf G S as follows. 



Qn 


B 


(id, • ■ • ,ji) 






{l,2,...,n- 


^1} 




(2,3,...,n + 


1) 


4(2) 


{1,2, 


. . , n, -n, . . 


,-2,-1} 


(2,3,.. 


.,n, -1,-n, . 


.,-3,-2) 


4(2) 
^2n 


{1,2,. 


. , n, —n, . . . 


-2,-1,0} 


(2,3,... 


,71, -1, -n,... 


,-3,-2,0) 


5(1) 


{1,2,. 


. , n, 0, —n, . 


..,-2,-1} 


(2,3,. 


. .,n,0,-n, .. 


,-3,-2) 


rW 


{1,2, 


..,n,-n,.. 


,-2,-1} 


(2,3,..., 


n, —1, — n, . . . 


-3,-2,-1) 




{1,2, 


..,n,-n,.. 


,-2,-1} 


(2,3, 


. . . , 71, -n, . . . 


-3,-2) 




{1,2,.. 


,n,0, -n, . . 


,-2,-1,0} 


(2,3,.. 


. , n, 0, — n, . . . 


-3,-2,0) 



Here should be understood as a symbol and not an empty set. By the definition d = n 
for aI^^ , d = 277 - 2 for dI^^ , d = 277 - 1 for A^^n- i,Bi^\d^ 2n for A^^n , and D^^l^ . 
We let 

W = {{... ,6„6,+i,...) e ••■ X B X S X ••• I 6, = 1 for |i| > 1} (1) 
be the set of states of the automaton at a fixed time. Setting 

J = {ju--- ,jd}, 

one finds 

'i3\{l,-l} ifsn^Bi^\Di'\Dl^l„ 



J 



B\{1} otherwise 
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For j G J (hence j ^ 1), we define 



ifjG {0,-1,0}, 
^ ^ ■ ifjGJ\ {0,-1,0}. 

The time evolution of the fl„-autoniaton T : W W takes the factorized form 

T^K,,.--K,,K,,. (3) 

To define the operator Kj-.W^W [j G J), one regards an element (. . . ,bi, h+i, . . . ) 
of W as the one dimensional array of boxes containing bi in the i-th box. (The (i + l)-th 
box to the right of the i-th box.) In what follows we regard the box containing 1 e B as 
an empty box. Thus, taking some b I) E B away from a box means the change of the 
local states b 1. Similarly, putting b{^ 1) d B into an empty box means the change 
1^6. Note that the boundary condition in (|l|) says that only finitely many boxes are 
non-empty. Under this convention the operator Kj is defined by the following algorithm 
which consists of 4 Steps: 

1. Replace each —1 by a pair within a box. 

2. Pick the leftmost j (if any) and move it to the nearest right box which is empty or 
containing just j* . (Boxes involving the pair j, j* are prohibited as the destination.) 

3. Repeat Step 2 for those j^s that are not yet moved until all of j's are moved once. 

4. Replace the pair within a box (if any) by —1. 



Remark 1. For g„ = An where no —1 is present in B and J, Step 1 and 4 become void 
and the above procedure shrinks to the well known ball-moving algorithm in the box-ball 
system Q. 

Remark 2. When j G { — 1,0,0}, one can have two j's within a box due to (||). As for 
those duplicated j's, it docs not matter which is left or right, but of course they should be 
distinguished according to whether they are moved or not yet moved. See K(ij in Example 

I 

Remark 3. A little inspection shows that K-i actually admits a simpler description which 
is essentially the ball-moving algorithm in the box-ball system: 

1'. Pick the leftmost —1 (if any) and move it to the nearest right empty box. 

2'. Apply Step 1' for those — I's that are not yet moved (if any). 

3'. Repeat Step 2' until all the —I's are moved once. 

Remark 4- Consider the g„ ^ A^P case and interpret the local states as follows: 



B 


Entry in the box 


1 


empty 




particle of color j 


-1 


bound state of j and j* 



Since the set 5 \ {1,-1} is invariant under the interchange j ^ j*, the state j* may 
be viewed as the anti-particle of j and vice versa. In this sense the bound state —1 is 
neutral, and especially and represent the 'self-neutral' particles that can still form 
a bound state with another one. Under this interpretation, the above algorithm for Kj 
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describes the motion of right-moving particles of color j seeking the empty box or a free 
partner, i.e., an anti-particle j* not yet paired with the other j's, to form a neutral bound 
state within a box. 

The algorithm for Kj can also be stated in terms of local rules, which we shall now 
explain. For j G J we introduce a map Lj : (Z>o) x B ^ B x (Z>o) as follows. Let 
b 

the diagram I denote Lj : i-> {b',l'). (j is attached to the horizontal line.) 

b' 



1. j ^ {1, 0, -1, 0} case. Assume I e Z>o, b e B\{j, -j, 1, -1}. 

3 -3 -j -1 1 1 



-1 -3 -j 



-l + l l + l- 



■/ 0- 



u 

l-^l 



-0 l—t^l 

b 



2. j £ {0, 0} case. Assume I G Z>o, b e B\{j, 1, -1}. 
-1 3 3 1 



l- 



^1 + 2 0- 



^1 i+i-^i+i i+2-^i 1—^0 0— 1^ 



3. j = —1 case. Assume / e Z>o, b £ B\{l,j}. 



3 1 
^l + l l + l-^l 
j 



/- 



u 



For 0„ = An'' we only have the third case with the j's understood as j £ {2,3,... ,n + l}. 



See also IHIK 



Given an automaton state (. . . ,bi, . . .) £ W, there exists an integer m such that 
bm' = 1 for all ra' < m owing to the boundary condition (|^). Fix any such m. Then the 
operator Kj : W ^ W maps (. . . , bi, bi+i, . . . ) to (. . . , c^, c^+i, . . . ), where Cm' — bm' for 
all m' < TO. The remaining c,„, Cm+i, . . . are determined by the composition of Lj's as 



bm bm+1 • • • 



(4) 



It is easy to see that the result is independent of the choice of to. The nonnegative integers 
I on the horizontal line stand for the number of color j particles on the carrier. The 
diagrams for Lj are viewed as the loading and unloading process of the color j particles 
when the carrier proceeds to the right. They match the algorithm stated before. Perhaps 
the above formulation of Kj using Lj is easier to program. 
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3 Examples 



Let us present examples of the factorized dynamics. We suppress the trivial left tail in 
the time evolution and only depict the part that corresponds to the composition of 
vertically. 

Example 1. 0„ = D'i\ T = K^K^K^K^^K^sK^^- 

-3 -2 1 -2 2 3 111 







■0- 



1- 



-0- 



■1- 



■0- 



■0- 



■0- 



-0- 



1 -1 



Ko 0-^1 



1-^0- 



111 

3 1-3 



-1111 



0-^1. 



0—^0- 



1 1 



KJ 



1-111 



K. 



1 

Kol — ^0—1^0 



1 1 



-0- 



-0- 



^0 — ^0 



-3 3 1 



-0 



1 



1 



-3-2 1 1 3-3 3 
T : (-3, -2, 1, -2, 2, 3, 1, 1, 1, . . .) ^ (1- -3, -2, 1, 1, 3, -3, 3, 1, • ■ ■)• 



Example 2. 0„ = B^^l T = K2K^KaK ^:,K ^2- 

-1 -3 3 -2 1 1 1 1 1 



K.l 



■1- 



-1- 



■1- 



■1- 



■0- 



-0- 



►0- 



-3 3 1 



-2111 







10-11 



kJ 0-^0 



0—^1- 



3—^1—^1 



1 1 1 



1—^0- 



1 1 1 



K^i — ^0 — ^0—^0 



3 3 3 

^0 — ^1—^1 



-2011 



1-^1- 



1 1 1 



-2 3 1 



K2\ 



■1- 



■0- 



-0- 



-0- 



-0- 



■0- 



■0- 



-0- 



►0- 



1 



2 1 1 -3 -2 -2 3 1 
T: (-1,-3,0,3,-2,1,1,1,1,1,...) ^ (1,2,1,1,-3,-2,-2,0,3,1,...). 
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Example 3. 0„ = Af\ T = K2K^iK^2Kq. 

-1 -2 2 -2 



1111 



kJ 0— F--2 



K-2 



2—^2- 



-1111 



1 1 



K-i 



1 1 



1 



-1 



Ko — ^0 



0—^0- 



-0- 



-0- 



■1- 



■1- 



-0- 



1 10 10-21-1-11 
T : (-1,-2,0,2,0,-2,1,1,1,1,...) ^ (1,1,0,1,0,-2,1,-1,-1,1,...). 
Along the Steps 1-4 in Section ||, the action of Kt/, in the above goes as follows. 



|-1 1 |-2 1 


1 ^ 1 


1 2 1 


1 ^ 1 


|-2 1 




|0 0| |-2 1 


1 1 


1 2 1 


1 ^ 1 


|-2 1 
















1 0| |-2 1 





1 2 1 





|-2 1 
















1 1 l~2 1 


|0 0| 


1 2 1 


|0 0| 


|-2 1 
















1 1 1 


1 ^1 


1 2 1 


|0 0| 


|-2 1 


U 1 














1 1 1 


1 ^1 


1 2 1 


1 ^1 


|-2 1 


U ^1 


1 1 1 


1 ^ 1 


1 ^ 1 


1 1 


|-2 1 


l-l 1 



Step 1 
Step 2 
Step 3 
Step 3 
Step 3 
Step 4 



As cautioned in Remark 2, one must distinguish the moved 0's and those not yet moved. 
Here we have marked the moved ones with underlines. 



4 Solitons 

To save the space we shall write, for example, 2^2(— 3)**-^ to signify the configuration of 

the local states 2,... ,2,-3,... ,—3 in a segment of an automaton state for nonnegative 
integers j/2 and ys, etc. In the sequel we assume yi, G Z>o for any b Cz B. For each g„ 
consider the following configurations and define v from them. 
A'n^ : (n-Kl)2'"+i ... 

En+l 
t=2 y»- 
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: (0)^»(-l)^-^ ("2)^- 

V = y9 + T,7=iiyt + y-^)- 
B^n^ : {-2)y-^ {-'d)y-'^ ... 

V = yo + YJ^=2iy^ + y-i)- 

Cl^'^ : {~l)y-^ (-2)s'-2 

« = YJl=iiyt + y-i)- 
oi^'' : {-2)y-^ {-3)y-^ 



. . . {~n)y-" ny 



{-n)y- 
{-n)y- 



Di^li ■■ (0)^''(-i)^-^ {-2)y-^ . 
v = yi!) + yo + Er=i(y» + y-^)■ 

In an element of W , the above configuration is called a soliton with amplitude v if it 
is surrounded by sufficiently many I's € B. (For an amplitude v soliton, 1" in its right 
suffices. ) The data {yt} is the internal label of a soliton. Under the time evolution T 
(H), the amplitude v solitons propagate stably to the rig ht with velocity v if they stay 
sufhciently away from configurations different from 1 G B\ HKTl | . It is a good exercise to 
check this claim by using the algorithm for T in Section |. Moreover, the follo wing fact s 
have been proved concerning the collisions of solitons with distinct amplitudes [ HKOTY | . 



' ... 2^= iy\ 2/0 e {0,1}, yi=y-i, 
...3y^2y\ yoe{0,l}, 

3y'2y\ yny-n = 0, 
Qvo „y„ ... 21/2 ivi^ y^^y^ e {0,1}, yi^y^i, 



1. The set of amplitudes remains invariant under the collisions. 

2. Two-soliton scattering rule is characterized by the combinatorial R of J7g(g„_i). 

3. Multi-soliton scattering factorizes into the two-soliton ones. 



Leaving the precise statements to [HKOTY], we here include an example of collision of 
two solitons. 

Example 4. g„ — C'^^K To make the space uniform, we denote —j by j for j ~ 1, 2, 3. 
Then the evolution of a two-soliton state under (0 < t < 5) is depicted as follows. 



t=0 
t=l 
t=2 
t=3 
t=4 
t=5 



. . 11122311111232111111111111111111111111 
. . 11111111223111232111111111111111111111 
. . 11111111111112311223211111111111111111 
. . 11111111111111111311122232111111111111 
. . 11111111111111111111311111222321111111 
. . 11111111111111111111111311111112223211 



One observes that the initial two solitons 12231 and 232 are scattered into the final 
solitons 131 and 22232 with a phase shift. In this way, the two-body scattering rule 
consists of the exchange of the internal labels of solitons and the phase shift. 



5 Relation to crystal theory 



The factorized dynamics in this paper is a translation of a result in | IIKT2 1 . By regarding 
the local states —1, . . . , —n here as 1, . . . ,n, B can be identified with the crystal Bi as a 
set. In eq.(33) of [ HKT2| , take the integer A; as = n — 1 for and k — n for the other 
g„'s. It leads to Ofc = 1 (cf. Table I, II and eq.(ll) therein) for all the 0„'s in question, 
which implies the boundary condition ... Ill ... . The time evolution in this case is given 
by T = Tk+d, where the data ik+d, ■ ■ ■ , ik+i in eq.(33) there reads 
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fin 


ik+d, • ■ • , ifc+1 






2, 3, . . . , n, n + 1 




.(2) o(l) 

4(2) ^(1) r)(2) 


0,2,3,. 
1,2,.. 


. , 71 — 1, n, 71 — 1, . . 

, 71 — 1, 71, 71 — 1, . . . 


.,3,2,0 
,2,1,0 


^(1) 


0,2,3, n- 


2, 71 — 1, 71, 71 — 2, 71 


-3, ...,3,2,0 



Under the convention 7^ = id, the operator Kj in this paper has emerged from the 
formula 

Kj^=Tk+r{rk+r-lT^ l<r<d. (5) 

We have verified that the right hand side yields the composition of ij,.'s described in 
Section ^ by an explicit calculation. Similar factorized dynamics can be formulated 
corresponding to the other boundary conditions than 1 e _B in (0) by choosing a different 
k in the above. However such variations do not affect the qualitative feature of the 
automata. 
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